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Description 

FIELD OF THE INVENTION 

[0001] This invention relates to process control sys- 
tems, and more particularly, to a system for determining 
whether a process parameter is in a steady or unsteady 
state through a use of the parameter's power spectral 
density. 

BACKGROUND OF THE INVENTION 

[0002] In controlling dynamic plant processes, it is of- 
ten necessary to know whether process variables are in 
a steady state or in an unsteady state. While, many proc- 
ess control systems monitor plant variables and com- 
pare them against predetermined set points, it is often 
more important to ioiow whether deviations of a plant 
variable, overtime, away from a setpoint are significant- 
ly different from normal. It is also important to be able to 
predict wrtiether a plant parameter output exhibits an in- 
cipient condition which may lead to an unsteady state. 
[0003] The prior art contains many teachings con- 
ceming plant parameter monitoring for process control 
applicmlons. In U.S. Patent 4,744,041 to Strunk et al., 
the steady state speed of a dc motor is detected through 
the use of fast Fourier transform analysis. Strunk et al 
employ a current sensor which measures the current in 
a test.motor and sends the sampled current signal as 
data to a computer The computer then samples and 
stores measured instantaneous current values at plural- 
ity of discrete times and performs a fast Fourier trans- 
form on the steady state current to determine its power 
spectral density. Based on the computed power spectral 
density, the speed of the motor is determined by detect- 
ing the frequency at which maximum power Is used. 
[0004] Kato et al., in U.S. Patent 4,303.979, discksse 
a system for monhoring frequency spectrum variations 
in output signals. The system initially determines an 
RMS average frequency value for an input signal. In ad- 
ditfon, it determines F\WiS values for each of a plurality 
of frequency subranges within the input signal. The sys- 
tem then monitors an Input test signal and obtains its 
RMS value and the average frequency of the overall in- 
put signal. If the determined reference and test frequen- 
cies differ substantially in their RMS and average fre- 
quency values, the RMS value and average frequency 
of an anomalous frequency component peak is cateu- 
lated. The average of the anomalous frequency compo- 
nent is then compared to boundary frequency values to 
determine in which frequency range the anomatous val- 
ue lies. The RMS value of the thus determined frequen- 
cy range and the average frequency are then employed 
to determine correctton parameters. 
[0005] Grassart. in U.S. Patent 4,965,757 discloses a 
process and device for decoding a received, encoded 
signal. The received signal is first filtered, sampled and 
digitized before being stored. Digitized samples of each 



successive signal block are transposed to the frequency 
domain by a fast Fourier transform. The thus computed 
spectra are compared with stored theoretical values for 
each possible code signal in order to identify the re- 
s ceived encoded signal. 

[0006] Schmidt, in U.S. Patent 3,883,726, discloses 
a fast Fourier transform algorithm computer that utilizes 
an attenuated input data window. An input buffer re- 
ceives input time samples, a cosine square attenuator 

10 superimposes a cosine squared shape on the input data 
samples, and then the resultant signal is passed to the 
fast Fourier analysis computer. A delay device, adders 
and an output buffer are provided for removing the effect 
of the attenuating input data window. 

IS [0007] Toda et al. in U.S. Patent 4,976.633, disclose 
a spectrum analyzer that displays spectrum data and 
power values of an RF or optk^al signal. An input signal 
is directed down one path where it is subjected to a 
spectrum analysis, and down a second path where its 

20 power value is determined. Display means indicates 
both the spectrum data and the power value. 
[0008] Murphy et al.. in U.S. Patent 5.087.873. dis- 
ctose apparatus for detecting a con'oslon state of buried 
metallic pipe. In Figs. 4-7. a spectrum analyzer receives 

2S input signals from a pair of magnetometers. The mag- 
netic field, as a function of frequency. Is determined by 
conducting a fast Fourier transform on the received sig- 
nals, with the resulting spectrum indfcating the ampli- 
tude and phase of the magnetk: field. Those values are 

30 used to determine the conditnn of the buried pipe. 
[0009] Demjanenko et al. In U.S. Patent 4.980,844 
teach a method for diagnosing mechanical conditions of 
a machine. The Demjanenko et al. procedure is compu- 
tationally complex and involves the creation of one or 

35 more sets of reference signatures (those reference sig- 
natures may be power spectrum data). Similar test data 
is then acquired and a comparison action occurs. The 
procedure averages the test and reference signatures, 
determines a distance therebetween (i.e. the Euclklean 

40 distance), compares that calculated distance against a 
computed threshold value, and classifies the distance 
as being normal or abnormal based upon the compari- 
son. The threshold is established as a function of the 
mean reference distance and reference and test stand- 

<s ard deviations. 

[0010] Toshiba Japanese patent application 
S62-245931 describes a procedure for determining ab- 
nornnalities in a rotating pump. A plurality of reference 
frequency spectra are stored in memory (called restric- 

so tion patterns), each reference frequency spectrum in- 
dicative of an operational state of the pump under dif- 
ferent load conditions. To conserve memory, a 'bias' 
procedure is employed to nxKlify a closest reference fre- 
quency spectrum to bring it into the load (I.e. power) 

55 range of a test frequency spectrum taken from the 
pump. If the test frequency spectrum exceeds limits de- 
termined from the 'biased' reference frequency spec- 
trum, an alarm is given. The Toshiba system thus re- 
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(c) storing said reference spectral data; 

(d) establishing current operational data for said 
process variable by sampling a current output sig- 

s nal for said process variable over a second period 
of time when said plant is in operation, each sam- 
pled current output signal being subjected to said 
first normalisation procedure to convert each sam- 
pled signal to a zero mean basis; 

10 

(e) processing the normalised sampled current out- 
put signal to derive current spectral data therefrom, 
said current spectral data manifesting an energy 
content of each of a plurality of frequency compo- 

is nents of the nomialised sampled cun^ent output sig- 
nal, each current spectral data value being subject- 
ed to said second normalisation procedure, wherein 
each spectral data value is divided by a sample time 
related value, a sample time being a window of time 

20 during which an output signal is sampled in step (d) 
in preparation for deriving current spectral data val- 
ues therefor; and 

(f) for each of a plurality of common frequency com- 
ponents of said reference spectral data and current 
spectral data, comparing nomnalised energy con- 
tents therefor and issuing a signal if said comparing 
indicates that the energy content of any frequency 
component of said current spectral data exceeds 

30 the energy content of a common frequency compo- 
nent of said reference spectral data by a predeter- 
mined energy threshold. 

DESCRIPTION OF THE DRAWINGS 

3S 

[0016] 

Fig. 1 is a block diagram of a system for performing 
the method of the invention; 

40 

Figs. 2a and 2b are high level flow diagrams illus- 
trating the method of the invention; 



quires storage of a plurality of reference frequency spec- 
tra and accurate modification of a chosen spectrum by 
application of a bias value thereto to enable its compar- 
ison with a test spectrum. Misadjustment of the bias can 
negate the effectiveness of the test system. 
[0011] Bartusiak et al. In PCT published Applk;ation 
WO 94/22025 describe a process signal detection tech- 
nique that is dependent upon power spectral densities. 
Reference data is accumulated while a plant operates 
at a steady state and Includes the energy content of 
each of a plurality of frecuency components. A current 
operational database is then established by sampling a 
process variable when the plant Is in operation. The 
sampled current output signal is analyzed to derive cur- 
rent data that includes the energy content of each of a 
plurality of frequency components. For each common 
frequency component of the reference data and current 
data, the system makes a comparison of the energy 
contents and issues a non-steady state signal if a differ- 
ence between the compared energy contents exceeds 
a predetermined limit. 

[0012] Accordingly, it is an object of this invention to 
provide a method for detemiining when a process vari- 
able output signal is in a steady state condition. 
[001 3] It Is another object of this invention to provide 
a method for analysis of a parameter output signal to 
indicate an incipient unsteady state. 
[0014] It is still another of this invention to provide a 
plant parameter output signal detection technique that 
. is dependent upon power spectral densities. 

. SUMMARY OF THE INVENTION 

[0015] A method for determining a current state of a 
process variable output signal from a plant and for de- 
temiining Whether the output signal is within acceptable 
limits, the method comprising the steps of: 

(a) establishing reference comparison data for saki 
process variable by sampling a reference output 
signal over a first period of time when said plant is 
operating at steady state, each sampled reference 
output signal being subjected to a first normalisation 
procedure to convert each sampled signal to a zero 
mean basis; 

(b) processing the normalised sampled output sig- 
nal to derive reference spectral data for said proc- 
ess variable, saki reference spectral data manifest- 
ing an energy content for each of a plurality of fre- 
quency components of the normalised sampled out- 
put signal, each reference spectral data value being 
subjected to a second normalisation procedure, 
wherein each spectral data value Is divided by a 
sample time related value, a sample time being a 
window of time during which an output signal is 
sampled in step (a) in preparatton for deriving ref- 
erence spectral data values therefor; 



Fig. 3 is a plot of a variation in flow over time in an 
exemplary plant; 

Fig. 4 is a semi-log plot of energy versus frequency 
derived from the signal shown in Fig. 3; 

Fig. 5 is a pk>t of flow vs. time showing a negative- 
going ramp disturbance; 

Fig. 6 is a semi-tog ptot of energy vs. frequency of 
the plot of Fig. 5 indicating the that the ramp distur- 
bance produces the most severe violation of plotted 
energy threshokis in the tow frequency range; 

Fig. 7 Is a plot of flow vs. tinne showing a "U' type 
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disturbance in the flow signal; 

Fig. 8 is a semi-log plot of energy vs. frequency of 
the plot of Fig. 7 indicating, at one energy threshold, 
that a disturbance exceeds the threshold. 

DETAILED DESCRIPTION OF THE INVENTION 

[001 7] The system and method for carrying out the in- 
vention determine whether a plant process is operating 
within a steady state or a non steady state condition. In 
brief, the system determines a reference power spectral 
density (PSD) of a process variable, during a time that 
the process variable is In a steady state. A power spec* 
tral density is a representation of an energy content of 
a signal as a function of its component frequencies and 
is computed via a fast Fourier transform (FFT) which 
converts a time-domain signal into its frequency-domain 
representation. 

(001 8] The signal from which the reference PSD is de- 
rived is subjected to a first normalization procedure to 
derive a zero mean basis for the signal. This normaliza- 
tion is accomplished by subtracting the mean process 
variable value of a time series of sample data values 
from each time series sample data value. A second nor- 
malization procedure is performed on the spectral data 
by dividing each power value in the spectrum by the 
number of time series signal samples used to compute 
the spectrum. The first normalization procedure enables 
a single stored reference PSD to be used as a reference 
in all subsequent test procedures and eliminates the 
need for plural reference spectra or for a biasing proce- 
dure to adjust the reference spectra The second nor- 
malization procedure enables the reference PSD to be 
compared to a current PSD, even in the event that the 
time windows differ during which the current and refer- 
ence PSD's are derived. 

[001 9] The reference PSD is compared with a current 
PSD derived when the process variable is in operation. 
The current PSD is also subjected to the first and second 
normalization procedures described above with respect 
to the reference PSD. If the current PSD is too energetic 
(either with respect to slow changes like ramps or with 
respect to fast changes like spikes) relative to the refer- 
ence PSD, the process variable is determined to be un- 
steady and such condition is signaled. 
[002(q Turning to Fig. 1 , flow monitors 1 0 and 1 2 con- 
tinuously monitor the state of flow in a pair of pipes 14 
and 16, respectively. Pipes 14 and 16 form portions of 
a plant whose process variables are continuously mon- 
itored to detemriine if any one has moved from a steady 
state to a unsteady state. Those skilled in the art will 
realize that the representation of Inputs from pipes 14 
and 16 is merely exemplary and that a plurality of other 
types of system variables (e.g. pressure, volume, tem- 
perature, etc.) can be nfK)nitored. The outputs from each 
of flow monitors 10 and 12 are fed to respectively con- 
nected anak>g to digital converters (A/D) 18, 20 whose 



outputs are, In turn, connected to a bus 22 that forms 
the main communicatkMi pathway in a control data 
processing system. 

[0021] A central processing unit (CPU) 24 is intercon- 
s nected with bus 22 and also is connected to A/D con- 
verters 18 and 20 to enable timed, sample signals to be 
derived therefrom. A read only memory (ROM) 25 is 
connected to bus 22 and contains a procedure for op- 
erating CPU 24 to monitor sensors 10. 12 and a proce- 

10 dure for enabling CPU 24 to perform a fast Fourier trans- 
form (FFT) of input data received from A/D converters 
1 8 and 20. A random access memory (RAM) 26 is con- 
nected to bus 22 and contains allocated memory for 
storing: raw input data from A/D converters 18 and 20; 

IS reference PSD data determined during a time that a 
process variable being monitored is In a steady state; 
and current PSD data that is detemiined when the proc- 
ess variable is being cun^ently monitored. 
[0022] The operation of the system of Fig. 1 will now 

20 be described in conjunction with the flow diagram shown 
in Figs. 2a and 2b. Initially, a steady state (i.e. reference) 
PSD is determined for a monitored process variable. 
The steady state PSD Is derived by initially time sam- 
pling a process variable output signal when the plant is 

2S operating in a steady state condition (box 30). To avoid 
contamination of the steady state PSD, the time sam- 
pled output signal is filtered to remove aperiodic signals 
therefrom (box 32). The filtered output samples are then 
normalized to a zero mean basis by subtracting the 

30 mean process variable value of a time series of sample 
data values from each time series sample data value. 
This eliminates any bias that is operating-state depend- 
ent (box 33). 

[0023] The filtered, normalized periodic sample sig- 
35 nals are then subjected to an FFT analysis by central 
processor 24, under control of a procedure read out from 
ROM 24 (box 34), A second normalization procedure 
eliminates any bias that results from the time window 
during which the time series samples were detected. 
4a The second nomnalization procedure divides each spec- 
tral frequency data value by the number of time samples 
so as to achieve a "tffne* normalizatkm of the frequency 
data. 

[0024] The result of the PSD analysis is a series of 
45 normalized energy values at frequencies (oo^ for the fil- 
tered steady state signal, each normalized energy value 
at a frequency {(ojj having an attribute associated there- 
with that is proportional to the energy contained in the 
particular frequency signal. Those frequencies and en- 
50 ergy value attributes are stored in RAM 26 as steady 
state (or reference) PSD data. 
[0025] At this point, a user-supplied multiplier factor 
is accessed that is to be applied to each energy value 
attribute of the steady state PSD. The multiplier factor 
ss enables the derivation of an energy threshold which is 
the amount that a process variable's cun'ent PSD must 
exceed a reference PSD to be considered in an un- 
steady state (box 36). 
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[0026] The system of Fig. 1 now switches to a ■cur- 
rent' monitor state wherein it monitors the outputs of 
sensors 10.1 2, etc. During the time that any one sensor 
is monitored, its output is converted by an FFT proce- 
dure to a PSO representation which is then compared s 
with the previousfy derived steady state PSD for the 
same sensor (i.e. is stored in RAM 26). CPU 24 com- 
mences the current monitor state by sampling output 
signals from a sensor (e.g. 10) during operation of the 
plant (Box 38). After sufficient samples have been ac- io 
cumulated, CPU 24 filters the sampled current output 
signal to remove aperiodic components (box 39), nor- 
malizes the filtered output signal to a zero-mean basis 
(identically as with the reference PSD) (box 40) and then 
computes a PSD(<0|) for the sampled current output sig- is 
nal using an FFT procedure (box 41 ). Each spectral fre- 
quency data value is now divided by the number of time 
series samples to achieve a 'time* normalization of the 
frequency data. Time nomnalizing of both reference and 
current PSD data enables those data to be compared 20 
even if different time windows are used during their de- 
tection. 

[0027] It Is now determined whether any cu rrent PSD 
((D|) is greater than its corresponding steady state PSD 
(coj), mult^iied by the user-inputted multiplier factor 
(boxes 42,44). If yes, the system outputs the value of g){ 
at which the maximum energy ratio is found and a signal 
is issued that the process variable is in a non steady 
state condition (box 46). 

[0028] If the decision indicated in decision box 44 is 30 
that no current PSD(co|) exceeds the multiplier factor 
times the steady state PSD(fioj). then It is detemnined that 
the sampled current output signal is in a steady state 
condition and no further action is required. At this point, 
the procedure terminates and CPU 24 may then monitor 55 
another current process variable signal and the proce- 
dure is repeated. 

[0029] Figs. 3-8 contain examples of the procedure 
described above. In Fig. 3. a process variable (flow) is 
plotted over a period of several hours (126 minutes). 40 
Each dot on the plot trace evidences a sampled flow val- 
ue that is input to the system. Prior to monitoring the 
current flow values, a reference flow during steady state 
of the same process variable was monitored and a ref- 
erence PSO was derived (e.g. dotted trace 100 in Fig. ^ 
4). Reference PSD trace 100 has been quantized into 
ten 'bins" so as to enable averaged comparisons to be 
made. The PSD for the signal of Fig. 3 is then calculated 
(trace 1 02) and is seen to exhibit considerable variations 
in power over the frequency spectrum. £0 
[0030] Turning to Fig. 5, a ramp disturbance 103 is 
plotted (flow vs. time). Furthermore, for exemplary pur- 
poses only, it is assumed that two multiplier factors (e. 
g. 1 .0 and 1 .5 ) have been utilized to derive energy PSD 
thresholds 108. 110 respectively. The conversion of ss 
ramp signal 103 to a PSD results in energy versus fre- 
quency trace 106 shown in Fig. 6. Note that the ramp 
disturbance 103 produces the most severe vblation of 
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both PSD energy thresholds 108 and 110 in the low fre- 
quency range. Under such condittons, an alarm is gen- 
erated. The system may be adjusted such that an alarm 
is inhibited from issuance if only energy threshold 108 
is exceeded, however the frequencies at which the en- 
ergy exceeds threshold 1 08 may be output for user mon- 
itoring of a possibly incipient instability nrxxle. 
[0031] In Fig. 7. a "U' disturbance 113 is depicted in 
a plot of flow versus time. Over the span of the time dur- 
ing which the flow is monitored in Fig. 7, the flow differ- 
ence Is not large and yet it is clear that a disturbance 
does occur. As shown by Fig. 8. such disturbance is 
picked up by comparison of the peak corresponding 
PSD trace 112 with energy threshold 114. 
[0032] It should be understood that the foregoing de- 
scription is only illustrative of the invention. \^rious al- 
ternatives and modifications can be devised by those 
skilled in the art without departing from the invention. 
Accordingly, the present invention is intended to em- 
brace all such alternatives, modifications and variances 
which fall within the scope of the appended claims. 

Claims 

1. A method for detemiining a current state of a proc- 
ess variable output signal from a plant and for de- 
termining whether the output signal is within accept- 
able limits, the method comprising the steps of: 

(a) establishing reference comparison data for 
saki process variable by sampling a reference > 
output signal over a first period of time when 
saw plant is operating at steady state, each 
sampled reference output signal being subject- 
ed to a first normalisation procedure to convert 
each sampled signal to a zero mean basis; 

(b) processing the normalised sampled output 
signal to derive reference spectral data for said 
process variable, said reference spectral data 
manifesting an energy content for each of a plu- 
rality of frequency components of the normal- 
ised sampled output signal, each reference 
spectral data value being subjected to a second 
normalisation procedure, wherein each spec- 
tral data value is divided by a sample time re- 
lated value, a sample time being a window of 
time during which an output signal is sampled 
in step (a) in preparation for deriving reference 
spectral data values therefor; 

(c) storing said reference spectral data; 

(d) establishing current operational data for 
said process variable by sampling a current 
output signal for said process variable over a 
second period of time when said plant is in op- 
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eration, each sampled current output signal be- 
ing subjected to said first nonmaiisation proce- 
dure to convert each sampled signal to a zero 
mean basis; 

5 

(a) processing the normalised sampled current 
output signal to derive current spectral data 
therefrom, said current spectral data manifest* 
ing an energy content of each of a plurality of 
frequency components of the normalised sam- io 
pled current output signal, each current spec- 
tral data vaiue being subjected to said second 
normalisation procedure, wherein each spec- 
tral data value is divided by a sample time re- 
lated value, a sample time being a window of . 
time during which an output signal is sampled 
in step (d) in preparation for deriving current 
spectral data values therefor; and 

(f) for each of a plurality of common frequency 20 
components of said reference spectral data 
and current spectral data, comparing nonmal- 
ised energy contents therefor and issuing a sig- 
nal if said comparing indicates that the energy 
content of any frequency component of said 2S 
current spectral data exceeds the energy con- 
tent of a common frequency component of said 
reference spectral data by a predetermined en- 
ergy threshold. 

30 

A method as recited In claim 1 wherein said sample 
time related value is a number of samples of an out- 
put signal talcen during said window of time. 

A method, as recited in claim 1 wherein step (f) fur- 3S 
ther indicates the common frequency for which said 
signal is issued. 

A method as recited in claim 1 wherein said analyz- 
ing steps (b) and (e) perfomi Fourier analyses of 40 
said sampled output signals to determine the fre- 
quency components and energy content of each 
thereof. 

The method as recited in claim 1 wherein said signal ^ 
Is only issued in step (f) if the energy content in a 
frequency component of said current data exceeds 
the energy content of a common frequency compo- 
nent of said reference data, multiplied by a set fac- 
tor, so 

The method as recited in claim 1 wherein, In step 
(f). said signal is issued after a ratio of the energy 
contents of each of said common frequency com- 
ponents Is determined and it is determined that at 
least one said ratio exceeds a predetermined factor. 



Patentanspruche 

1. Verfahren zum Bestimmen des aktuellen Zustands 
eines Ausgangssignals einer Prozessvariablen aus 
einer Aniage und zum Bestimmen, ob das Aus- 
gangssignal innerhaib alaeptabler Grenzen liegt, 
bel dem in Stufen: 

(a) f Or die Prozessvariable Referenzvergleichs- 
daten emnittelt werden. indem Dber einen er- 
sten Zeitraum. wenn die Aniage im stationaren 
Zustand arbeitet, ein Referenzausgangssignal 
abgetastet wird, wobei jedes abgetastete Re- 
ferenzausgangssignal einem ersten Normaii- 
sierungsverfahren unten/vprfen wird, um jedes 
abgetastete Signal in eine Null-Mittelwertbasis 
umzuwandein, 

(b) das normalisierte. abgetastete Ausgangssi- 
gnal zur Ableitung von Referenzspektrafdaten 
fur die Prozessvariable verarbertet wird, wobei 
die Referenzspektraldaten fur jede einer Viel- 
zahl von Frequenzkomponenten des normali- 
sierten, abgetasteten Ausgangssignals einen 
Energiegehalt anzeigen, wobei jeder Referenz- 
spektraldatenwert einem zweiten Normalisie< 
njngsverfahren unterworfen wird, wobei Jeder 
Speklraldatenwert durch einen mit der Abtast- 
zeit verknDpften Wert geteilt wird, wobei Abt- 
astzelt ein Zeitfenster ist. wahrenddessen ein 
Ausgangssignal in Stufe (a) bei der \A3rt}eref- 
tung zur Ableitung von Ref erenzspektrakJaten- 
werten f Or dieses abgetastet wird, 

(c) die ReferenzspektraMaten gespeichert wer- 
den. 

(d) aktuelle Betrlebsdaten fur die Prozessvaria- 
ble ermittelt werden, indem ein aktuelles Aus- 
gangssgignal fOr die Prozessvariable Qber ei- 
nen zweiten Zeitraum, wenn die Aniage in Be- 
trieb ist, abgetastet wird, wobei jedes abgeta- 
stete aktuelle Ausgangssignal dem ersten Nor- 
malisieningsverfahren unterworfen wrd, um je- 
des abgetastete Signal in eine Null-Mittelwert- 
basis umzuwandein, 

(e) das normalisierte. abgetastete aktuelle 
Ausgangssignal zur Ableitung aktueller Spek- 
trakjaten von diesem verarbeitet wird. wobei 
die aktuellen Spektratdaten etnen Energiege- 
halt von jeder von einer Vielzahl von Frequenz- 
komponenten des nonmalisierten. abgetaste- 
ten aktuellen Ausgangssignals anzeigen. wo- 
bei jeder aktuelle Spektrakiatenwert dem zwei- 
ten Normalisierungsverfahren unterworfen 
wird, wobei jeder Spektrakiatenwert durch ei- 
nen mit der Abtastzeit verknQpften Wert geteilt 
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wird. wobei Abtastzeit ein Zeitfenster ist. wah- 
renddessen eIn Ausgangssignal in Stufe (d) bei 
der Vorbereitung der Ableitung von aktuelien 
Spektraldatenwerten fur dieses abgetastet 
wird, und s 

(f) fur jede einer Vielzahl von gemeinsamen 
Frequenzkomponenten der Referenzspektral- 
daten und der aktuelien Spektrakiaten die nor- 
malisierten Energiegehalte fur diese verglichen io 
werden und ein Signal ausgegeben wird, falls 
der Vergleich zeigt, dass der Energiegehalt ir- 
gendeiner Frequenzkomponente der aktuelien 
Spektratdaten den Energiegehalt einer ge- 
meinsamen Frequenzkomponente der Refe- is 
renzspektraldaten um einen vorbestimmten 
Energieschweilenwert uberschreitet. 

2. Verfahren nach Anspruch 1, bei dem der mit der 
Abtastzeit verknOpfte Wert eine Reihe von Abtast- 20 
werten eines Ausgangssignals ist. die wahrend des 
Zeitfensters aufgenommen werden. 

3. Verfahren nach Anspruch 1 , bei dem Stufe (0 f emer 

die gemetnsame Frequenz anzeigt, fur die das Si- 2S 
gnal ausgegeben wird. 

4. Verfahren nach Anspruch 1 , bei dem die analysie- 
renden Stufen (b) und (e) Fourier-Analysen der ab- 
getasteten Ausgangssignale durchfuhren. um die 30 
Frequenzkomponenten und den Energiegehalt von 
jedem derselben zu bestimmen. 

5. Verfahren nach Anspruch 1. wobei das Signal in 
Stufe (f) nur ausgegeben wird, falls der Energiege- 3S 
halt in einer Frequenzkomponente der aktuelien 
Daten den Energiegehalt einer gemeinsamen Fre- 
quenzkomponente der Referenzdaten. multipliziert 

mit einem vorgegebenen Faktor. uberschreitet. 

40 

6. Verfahren nach Anspruch 1, wobei in Stufe (f) das 
Signal ausgegeben wird, nachdem ein Verhaltnis 
der Energiegehalte von jeder der gemeinsamen 
Frequenzkomponenten bestimmt worden ist und es 
bestimmt worden ist, dass zumindest ein solches 45 
Verhaltnis eInen vorgegebenen Faktor Qberschrei- 

tet 



Revendicatlons so 

1. Proc^dd pour determiner un dtat en cours d'un si- 
gnal de sortie de variable de processus d'une ins- 
tallatton et pour determiner si le signal de sortie ren> 
tre k I'interieur de limites acceptables. le procddd ss 
comprenant les dtapes suivantes : 

(a) retablissement de donndes de comparai- 



son de rdf^rence pour ladite variable de pro- 
cessus en echantillonnant un signal de sortie 
de reference au cours d'une premiere pdriode 
de temps pendant laquelle ladite installation 
fonctionne k I'etat stable, chaque signal de sor- 
tie de reference echantillonnd etant soumis k 
une premiere procedure de normalisation pour 
convertir chaque signal echantillonn6 sur une 
base de moyenne zdro ; 

(b) le traitement du signal de sortie echantillon- 
r\6 normalise pour deriver des donndes spec- 
trales de reference pour ladite variable de pro- 
cessus, lesdites donnees spectrales de refe- 
rence manifestant un contenu d'energie pour 
chacune d'une pluralite de composantes de fre- 
quence du signal de sortie echantillonne nor- 
malise, chaque valeur de donnees spectrales 
de reference etant soumise ^ une deuxieme 
procedure de normalisation, dans laquelle cha- 
que valeur de donnees spectrales est divisee 
par une valeur liee k un temps d'echantiltonna- 
ge. un temps d'echantilionnage etant unefenS- 
tre de temps pendant laquelle un signal de sor- 
tie est echantillonne au cours de fetape (a) en 
vue d'une preparatbn pour deriver des valeurs 
de donnees spectrales de reference pour celui- 
ci ; 

(c) la memorisation desdites donnees spectra- 
les de reference ; 

(d) retablissement de donnees de fonctbnne- 
ment en cours pour ladite variable de proces- 
sus en echantiltonnant un signal de sortie en 
cours pour ladite variable de processus au 
cours d*une deuxieme periode de temps lors- 
que ladite installation est en fonctionnement, 
chaque signal de sortie en cours echantillonne 
etant soumis k ladite premiere procedure de 
normalisation pour convertir chaque signal 
echantillonne sur une base de nrK>yenne zero : 

(e) le traitement du signal de sortie en cours 
echantillonne normalise afin de deriver des 
donnees spectrales en cours k partir de celui- 
ci, lesdites donnees spectrales en cours mani- 
festant un contenu d'energie de chacuno d'une 
pluralite de composantes de frequence du si- 
gnat de sortie en cours echantillonne nornnali- 
se, chaque valeur de donnees spectrales en 
cours etant soumise k ladite deuxieme proce- 
dure de normalisation, dans laquelle chaque 
valeur de donnees spectrales est divisee par 
une valeur liee k un temps d'echantilionnage, 
un temps d'echantilionnage etant une fendUe 
de temps pendant laquelle un signal de sortie 
est echantilbnne au cours de I'etape (d) brs de 
la preparatton pour deriver des valeurs de don- 
nees spectrales en cours pour celui-ci ; et 

(f) pour chacune d'une pluralite de composan- 
tes de frequence communes desdites donnees 
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spectrales de r^f^rence et desdites donn^es 
spectrales en cours, la comparaison de conte- 
nus d*dnergie normalises pour ce!les-ci et la dd- 
livrance d*un signal si ladite comparaison indi- 
que que le contenu d'energte d'une quelconque s 
composante de fr^uence desdites donn^es 
spectrales en cours ddpasse le contenu d'dner- 
gie d'une composante de frequence commune 
desdites donndes spectrales de reference d'un 
seuil d*6nergie prMdtermind. to 

2. Proc^dd selon la revendication 1 . dans lequel ladite 
valeur lide d un temps d'dcliantillonnage est un 
nombre d'^cliantillons d'un signal de sortie extrait 
durant ladite fenStre de temps. 

3. Proc^dd selon la revendication 1 , dans lequel I'^ta- 
pe (0 indique de plus la frequence commune pour 
laquelle ledit signal est d^livrd. 

20 

4. Proc6d6 selon la revendication 1 , dans lequel les- 
dites 6tapes d'anaiyse (b) et (e) effectuent des ana- 
lyses de Fourier desdits signaux de sortie tehan- 
tillonnds alin de determiner les composantes de fre- 
quence et le contenu d'dnergie de chacune de eel- ^5 
les-ci. 

5. Procede selon la revendication 1 , dans lequel ledit 
signal est deiivre uniquement lors de retape (f) si le 
contenu d'energie dans une composante de frd- 30 
quence desdites donnees en cours depasse le con- 
tenu d*energ!e dune composante de frequence 
commune desdites donnees de reference, muitiplie 
par un facteur determine. 

3S 

6. Precede selon la revendication 1 , dans lequel, au 
cours de retape (f), ledit signal est deiivre aprds 
qu'un rapport des contenus d'energie de chacune 
desdites composantes de frequence commune ait 
ete detemnine, et il est detemnine qu'au moins fun 40 
desdits rapports depasse un facteur predetermine. 
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FIG.2A 



STEADY- STATE PSD 
DETERMINATION 



SAMPLE PROCESS VARIABLE 
OUTPUT SIGNAL WHEN PLANT 
OPERATING AT STEADY STATE 



30 



FILTER OUTPUT SIGNAL TO 
REMOVE APERIODIC COMPONENTS 
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NORMALIZE OUTPUT SIGNAL TO 
ZERO MEAN BASIS . (SUBTRACT 
MEAN VALUE FROM EACH TIME 
SERIES DATA) 
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PSD(nj) FOR OUTPUT 



•COMPUTE 
SIGNAL USING FFT 
•NORMALIZE PSD(nj) BY 
NUMBER OF SAMPLES', 

-LABEL AS REFERENCE PSD 
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ACCESS MULTIPLIER FACTOR 36 
TO FIG.2B 
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F1G.2B 



CURRENT PSD 
DETERMINATION 



SAMPLE PROCESS VARIABLE 
OUTPUT SIGNAL DURING 
CURRENT PLANT OPERATION 
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FILTER OUTPUT SIGNAL TO 
REMOVE APERIODIC COMPONENTS 



NORMALIZE OUTPUT SIGNAL TO 
ZERO MEAN BASIS .(SUBTRACT 
MEAN VALUE FROM EACH TIME 
SERIES DATA) 



•COMPUTE PSD (rj. ;} FOR OUTPUT 
SIGNAL USING FFT; 
■NORMALIZE PSD(rii)BY 
NUMBER OF SAMPLES; 
■LABEL AS CURRENT PSD. 
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FOR EACH COMMON FREQUENCY 

IN PSD(ni) AND PSD(iij), FIND RATIO 

PSD u-Vi)/ MULTIPLIER FACTOR- PSD (/ij) 
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YES 



IS ANYPSDlAj T 
> 

^ULT. FACTOR • PSD( A j ) 
? 



STEADY STATE 



OUTPUT: 



•A-, AT WHICH MAXIMUM 
RATIO OCCURS 

-ALARM SIGNAL 
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